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Surface  destruction  and  change  of  solids  properties 
under  the  plasma  influence. 

A.L.  Bondareva,  G.I.  Zmievskaya 

Keldysh  Institute  of  Applied  Mathematics  Russian  Academy  of  Sciences 

Interactions  between  ions  beams  and  solids  or  plasma  and  solids  result  in  changing  of  solids 
properties  as  well  as  ions  beams  or  plasma  characteristics.  The  influence  of  blistering  on 
solids  properties  are  examined  in  the  paper.  The  stochastic  model  of  blistering  is  presented. 
The  evolution  of  distribution  functions  of  bubbles  versus  bubbles  sizes  and  coordinates  in 
lattice  is  received.  The  depth  dependences  of  porosity  and  tensions,  which  are  created  by 
bubbles,  are  calculated  using  the  distribution  functions. 


1.  Introduction. 

The  study  of  phenomena  of  interaction  between  plasma 
or  ions  beams  and  solids  surface  is  very  important  for 
technologies.  It  concerns  with  cosmophysics,  nanotech- 
nology and  controlled  thermonuclear  reactor.  The  phe- 
nomena result  in  change  of  surface  relief  and  properties 
of  materials.  Experiments  in  space  and  thermonuclear 
experiments  are  very  difficult  and  expensive.  In  this 
connection,  computer  simulation  of  such  processes  is 
very  important. 

This  paper  is  devoted  to  computer  simulation  of  fluctu- 
ation stage  of  high-temperature  blistering.  Blistering 
is  phenomenon  of  development  of  gas  bubbles  in  solids 
surface  layer.  Such  gas  bubbles  are  named  blisters. 
Blistering  results  from  irradiation  of  solids  surface  by 
ions  of  bad  soluble  gases.  The  fluctuation  stage  is  very 
short,  but  it  determines  further  development  of  blis- 
tering. The  duration  of  the  stage  is  10-4  sec  approx- 
imately. Bubbles  sizes  are  approximately  10A  during 
this  stage. 

The  behaviour  of  gas-vacancy  pores  in  Ni  crystal 
lattice  under  the  influence  of  He  ions  with  energy  from 
1 keV  to  35  keV  and  radiation  dose  from  1016  to  1019 
ions/cm2  is  examined  in  this  paper.  Such  terms  can 
arise  in  solar  wind  and  thermonuclear  reactor.  The 
temperature  of  material  T is  0.4-Tmen  < T < 0.6  Tmeu 
, Tmeit  is  melting  temperature  of  Ni. 

2.  Stochastic  model  of  blistering  and 
equations. 


solution  of  stochastic  differential  equations.  Scheme 
of  splitting  on  physical  processes  and  coordinates  was 
used  for  solution  of  examined  problem.  The  system  of 
four  stochastic  differential  equations  was  received  af- 
ter the  splitting.  Increase  of  size,  migration  on  lattice, 
fusion  of  bubbles,  exit  of  blisters  on  surface  and  bub- 
bles destruction  on  it  are  examined  in  this  paper.  The 
stochastic  equation  for  size  change  is  following 
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g is  bubble  size,  A<k(g,f,t)  is  Gibbs  potential, 
D(g,t)  is  diffusion  coefficient  in  sizes  space,  Wg{t)  is 
Wiener  process,  T*,  is  duration  of  fluctuation  stage. 

= 0,  gc r is  ’’critical”  size  of  blister. 

A$(gOmm!  T ) f)  = A$(gomaxi  f,  t)  — A$(gcr.,  r ,t)  — kT , 

90 min  < 9cr  < 90max- 

A $(g,  f,  t)  and  D(g,  t)  nonlinearly  depend  from  bubble 
size.  Heterogenous  condensation  on  bubbles  surfaces, 
difference  between  chemical  potential  of  phases,  surface 
tension  on  bubbles  surface,  elastic  reaction  of  solids 
lattice,  nonequivalence  of  bubbles  position  in  lattice, 
break  of  connections  of  crystal  lattice  are  considered 
using  A $(g,r,t)  and  D(g,t). 

The  stochastic  equation  for  depth  is  following 


The  stochastic  model  of  fluctuation  stage  of  high- 
temperature  blistering  had  been  put  forward  by  au- 
thors [1,2].  Brownian  motion  model  of  particles  with 
variable  mass  was  accepted  as  a basis.  Blistering  inves- 
tigation deals  with  computer  simulation  of  superposi- 
tion of  the  stochastic  processes  of  bubbles  sizes  change 
and  bubbles  migration  in  crystalline  lattice.  It  is  pos- 
sible because  processes  of  increase  of  bubbles  sizes  and 
bubbles  migration  in  lattice  have  appreciably  different 
time  scales.  Solution  of  kinetic  equations  for  Brow- 
nian motion  model  is  non-linear  problem.  Therefore 
authors  use  stochastic  analog  method  to  solve  these  ki- 
netic equations.  Main  idea  of  this  method  is  change 
of  kinetic  equations  by  their  stochastic  analogs  and 


z(t)=z(t0)+  f Hz(T,x(r),y(T),z{T))dT+ 
Jt  o 

+ f <r(T,z(r))dWz{T) 


Hz{t,x{t),v{t),z{t)) 
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1 dU  1 8DZ 

7 Mg  dz  2 dz 

U = Ubb  + Ubs 


Wz  is  Wiener  process,  z is  depth  of  bubble  from  surface 
under  irradiation,  Mg  is  blister  mass,  x,  y,  z are  bubble 
coordinate  in  crystal  lattice,  Dz  is  diffusion  coefficient, 
Ubb  is  interaction  potential  between  blisters  by  phonons 
and  oscillation  of  electron  density,  Ubs  is  potential  of 
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interaction  between  bubbles  and  surface. 

0.6  - br  + aT  cos(cr(fi  — fj)) 

\ri~fj\3 

(xj  - Xj)4  + (yj  - yjY  + (zj  - zj)4 
\n-?iV  J 

ar,  br  Cr  are  model  coefficients,  f,  and  fj  are  radius- 
vector  of  two  interactive  bubbles.  Two  bubbles  can 
fusion  if  fusion  condition  is  carried  out.  The  fusion 
condition  is  \n  -fj | < + p]/3)  + A/.  A/  is 

model  parameter.  0 < A / < a,  a is  lattice  parameter, 
rjie  is  radius  of  He  atom. 

Authors  modified  Artem’ev  method  to  solve  the 
stochastic  equations  and  used  new  method  to  solve  sys- 
tem of  four  concerned  with  each  other  stochastic  equa- 
tions with  nonlinear  coefficients. 

The  developed  model  allows  to  receiv  the  following 
characteristics:  distribution  function  of  bubbles  from 
bubbles  size  and  position  on  crystal  lattice,  mathe- 
matics] expectations  and  dispersions  of  bubble  size 
and  distance  from  surface,  evaluations  of  porosity  of 
solids  layers  and  tensions  in  layers. 

3.  Conclusions. 

1.  According  to  classical  theory  bubbles  destroy  if 
their  sizes  are  less  than  gcr  and  bubbles  increase 
if  their  sizes  are  more  than  gcr.  However,  numer- 
ical experiments  demonstrate  that  blisters  can 
increase  if  their  sizes  are  more  than  gamin  and 
less  than  gcr.  Bubbles  can  destroy  if  their  sizes 
belong  to  interval  from  gcr  to  gomax-  The  part 
of  increased  blisters  with  sizes  € [<7omini  gcr]  is 
larger  than  the  part  of  destroyed  blister  with  sizes 

€ [<?cr!  ffOmnz]. 

2.  The  bubbles  migration  into  the  surface  under  ir- 
radiation if  bubble  radius  is  less  than  5 A,  other- 
vase  bubbles  stop; 

3.  Layer-like  structure  of  near-surface  layer  is  dis- 
covered. Layers  with  accumulation  of  bubbles  al- 
ternate with  layers  where  number  of  blisters  is 
small; 

4.  The  greatest  porosity  is  observed  on  depths  of 
~ 0.85  • Rj,  and  ~ 0.35  • Rp,  Rp  is  middle  depth 
of  projection  run;  The  porosity  of  solid  layers 
from  depth  at  time  of  finish  of  fluctuation  stage 
is  shown  on  Figure  1.  The  depths  of  two  biggest 
maximums  of  porosity  correspond  to  depths  of 
development  of  two  kinds  of  blisters.  The  depth 
~ 0.85  ■ Rp  conforms  to  depth  of  ’’large"  bubbles. 
And  the  depth  ~ 0.35  • Rp  agrees  with  "small” 
bubbles. 

5.  Rate  of  increasing  of  bubble  size  reduces  if  blister 
size  is  bigger  than  12  A.  This  phenomenon  con- 
nects with  great  number  of  connections  breaks  in 
crystal  lattice; 
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Figure  1:  The  porosity  of  solid  layers  from  depth  at 
time  of  finish  of  fluctuation  stage  is  represented  on  this 
picture.  The  depth  from  surface  (z)  is  measured  in 
lattice  parameters.  0 is  surface.  400  corresponds  to 
2 ■ Rp.  The  porosity  is  considered  as  ratio  of  layer 
porosity  from  porosity  of  all  sample  at  initial  time. 


Figure  2:  This  figure  shows  the  tensions  in  solid  layers 
at  finish  time.  The  tension  is  measured  in  2 • 109 Pa. 
The  depth  from  surface  (z)  is  measured  in  lattice  pa- 
rameters. 0 is  surface.  400  corresponds  to  2 • Rp. 

6.  Tensions  in  solid  layers  at  finish  time  are  shown 
on  Figure  2.  Largest  tensions  are  observed  on 
depths  ~ 0.85  • Rp  and  ~ 0.35  • Rp.  Stress  is 
pronounced  near  surface. 

7.  Distribution  functions  of  bubbles  from  sizes  and 
coordinates  in  lattice  arc  nonequilibrium. 
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